Propagation of 1D Waves in Regular Discrete Heterogeneous Media: A Wigner Measure Approach

In this article, we describe the propagation properties of the one-dimensional wave and transport equations with variable coefficients semi-discretized in space by finite difference schemes on non-uniform meshes obtained as diffeomorphic transformations of uniform ones. In particular, we introduce and give a rigorous meaning to notions like the principal symbol of the discrete wave operator and the corresponding bi-characteristic rays. The main mathematical tool we employ is the discrete Wigner transform, which, in the limit as the mesh size parameter tends to zero, yields the so-called Wigner (semi-classical) measure. This measure provides the dynamics of the bi-characteristic rays, i.e., the solutions of the Hamiltonian system describing the propagation, in both physical and Fourier spaces, of the energy of the solution to the wave equation. We show that, due to dispersion phenomena, the high-frequency numerical dynamics does not coincide with the continuous one. Our analysis holds for C0,1(R)-coefficients and non-uniform grids obtained by means of C1,1(R)-diffeomorphic transformations of a uniform one. We also present several numerical simulations that confirm the predicted paths of the space-time projections of the bi-characteristic rays. Based on the theoretical analysis and simulations, we describe some of the pathological phenomena that these rays might exhibit as, for example, their reflection before touching the boundary of the space domain. This leads, in particular, to the failure of the classical properties of boundary observability of continuous waves, arising in control and inverse problems theory.